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Abstract 


The  intent  of  this  research  project  is  to  explore  alternatives  to  conventional  dipole  inversion  for 
extracting  target  features  from  multi-axis  EMI  sensor  array  data.  Conventional  dipole  inversion 
searches  for  the  target  parameter  values  (location  and  polarizabilities)  which  minimize  the 
difference  between  measured  signals  and  those  calculated  using  the  dipole  response  model.  Here, 
we  consider  an  alternative  approach  that  seeks  to  determine  the  parameter  values  which 
minimize  an  objective  function  based  on  the  dispersion  in  estimates  of  the  target’s  polarizability 
using  different  combinations  of  transmitters  and  receivers.  A  fortiori  the  approach  allows  for 
direct  calculation  of  the  uncertainty  in  the  estimated  polarizability. 

This  interim  report  documents  results  on  the  convergence  properties  of  a  downhill  simplex  based 
algorithm  for  detennining  a  target’s  location  (and  polarizabilities)  using  this  approach. 
Specifically  we  examined  convergence  of  the  algorithm  for  data  collected  with  the  2x2  TEM 
array  at  the  fonner  Camp  Beale  in  201 1  and  found  no  impact  due  to  signal -to-noise  ratio  (SNR) 
and  background  leveling  effects.  However,  the  minimum  polarizability  dispersion  does  vary 
systematically  with  SNR  -  targets  with  higher  SNR  tend  to  have  less  uncertainty  (dispersion)  in 
the  polarizability  estimate  than  those  with  lower  SNR. 

Testing  shows  that  polarizability  dispersion  based  inversion  can  produce  more  accurate 
polarizabilities  than  conventional  dipole  inversion  in  some  cases.  However,  the  fraction  of 
targets  in  the  ESTCP  classification  demonstrations  that  cannot  be  accurately  analyzed  using 
conventional  dipole  inversion  may  be  small  enough  (<5%)  that  modest  improvements  in 
calculating  the  polarizability  are  likely  to  have  very  little  effect  on  classification  perfonnance. 
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Objective 


The  Strategic  Environmental  Research  and  Development  Program  (SERDP)  and  the 
Environmental  Security  Technology  Certification  Program  (ESTCP)  have  invested  heavily  in 
developing  advanced  electromagnetic  induction  (EMI)  sensor  technology  capable  of  providing 
the  infonnation  needed  to  reliably  distinguish  between  buried  munitions  items  and  metallic 
clutter  in  the  ground  such  as  munitions  fragments,  scraps  of  exploded  metal  objects,  cultural  or 
agricultural  artifacts  and  debris,  etc.  The  processing  and  analysis  techniques  that  are  used  to 
extract  the  information  needed  to  reliably  classify  targets  as  munitions  or  clutter  from  data 
collected  with  these  sensors  were  originally  developed  for  use  with  an  older,  less  sophisticated 
generation  of  EMI  sensors  such  as  the  Geonics  EM61.  The  newer  sensors  provide  a  richer,  multi¬ 
dimensional  view  of  the  target  which  is  not  fully  exploited  by  the  processing  schemes  that  are 
currently  used. 

The  objective  of  this  project  is  to 

a)  Improve  understanding  of  the  relationships  between  multi-axis  EMI  sensor  array  data  and 
the  intrinsic  features  of  target  response 

b)  Develop  new  processing  approaches  for  multi-axis  EMI  sensor  data  that  make  full  use  of 
the  capabilities  available  with  newly  developed  EMI  sensor  technology,  and 

c)  Improve  procedures  for  identifying  targets  that  cannot  be  reliably  classified. 

Our  basic  approach  involves  exploring  new  processing  techniques  that  can  be  used  to  calculate 
the  standard  target  features  used  for  classification  (principal  axis  polarizabilities),  but  that  do  not 
rely  on  conventional  dipole  inversion.  In  conventional  dipole  inversion,  given  a  set  of 
measurements  of  the  EMI  response  we  simultaneously  search  out  the  location,  orientation,  and 
principal  axis  polarizabilities  that  produce  signals  that  best  match  the  measured  target  response. 
With  the  new  multi-axis  technology,  the  transmit  (Tx)  and  receive  (Rx)  fields  can  be  combined 
in  simple  ways  to  selectively  excite  and  observe  the  target's  response  in  a  set  of  orthogonal 
directions.  This  can  be  accomplished  using  appropriate  linear  combinations  of  the  data  recorded 
in  the  sensor's  various  transmit  and  receive  channels,  and  allows  the  polarizabilities  to  be 
calculated  directly  once  the  target's  location  has  been  detennined.  Different  combinations  of 
transmit  and  receive  coils  correspond  to  different  views  of  the  target,  and  we  need  only  find  the 
target  location  for  which  those  different  views  produce  consistent  results  for  the  target 
parameters. 

This  interim  report  documents  results  on  the  convergence  properties  of  our  algorithm  for 
detennining  a  target’s  location  (and  polarizabilities)  using  this  approach.  Specifically  we 
examine  the  effects  of  signal-to-noise  ratio  (SNR)  and  background  leveling  on  convergence 
using  data  collected  at  the  fonner  Camp  Beale  in  201 1. 


1 


Technical  Approach 


Electromagnetic  induction  is  the  process  whereby  changing  magnetic  fields  create  currents  in 
electrically  conducting  objects.  Transient  electromagnetic  (TEM)  sensors  create  a  magnetic  field 
by  running  a  current  through  a  coil.  When  the  current  is  shut  off  this  primary  field  shuts  down, 
inducing  currents  in  any  nearby  metallic  objects.  The  induced  currents  quickly  decay  as  they 
diffuse  through  an  object.  TEM  sensors  measure  the  voltage  induced  in  a  coil  by  the  decaying 
magnetic  field  associated  with  these  eddy  currents. 

Dipole  Response  Model 

In  the  standard  dipole  response  model,  we  represent  the  target  response  by  an  induced  dipole 
moment  m(t)  which  is  proportional  to  the  primary  field  Ho  at  the  target  location.  The 
proportionality  factor  is  the  target’s  magnetic  polarizability  tensor  B(t).  The  elements  Bij  of  B 
form  a  3x3  array  of  dipole  moment  components  in  the  i  coordinate  directions  induced  by 
excitation  in  the  j  coordinate  directions  such  that 

mx(t)  =  Bxx(t)H0x  +  Bxy(t)H0y  +  Bxz(t)H0z, 

my(t)  =  ByX(t)H0x  +  Byy(t)H0y  T  ByZ(t)HQz, 

mz  00  =  Bzx(t)H0x  +  Bzy(t)H0y  +  Bzz(t)H0z. 

Receive  coils  measure  the  signal  from  the  induced  dipole  field.  The  sensitivity  to  the  induced 
dipole  components  mimics  the  field  that  would  be  produced  by  currents  flowing  through  the 
receive  coil.  Both  the  primary  field  and  the  receive  coil  response  are  calculated  as  Biot-Savart 
integrals  around  the  coils.  Writing  T  for  the  field  of  the  transmit  coil  and  R  for  the  pseudofield  of 
the  receive  coil  (including  receiver  gain  factors),  the  signal  from  the  transmit/receive  coil  pair  is 
then  simply 


S  =  R  •  BT. 

B  is  a  second  rank  tensor  that  depends  on  the  size,  shape  and  material  properties  of  the  target,  as 
well  as  its  orientation  relative  to  the  x,  y,  z  coordinate  directions.  The  target  orientation  relative 
to  the  x,  y,  z  coordinate  axes  is  specified  by  a  set  of  roll,  pitch  and  yaw  angles  \|/,  9  and  (j).  If  the 
coordinate  system  is  rotated  by  the  angles  \\i,  9  and  (j)  into  alignment  with  the  principal  axes  of 
the  target,  then  the  off  diagonal  elements  of  B  go  to  to  zero  (By  =  9  for  i  ^  j),  and  the  diagonal 
elements  correspond  to  the  eigenvalues  pi  of  B.  These  are  referred  to  as  the  target’s  principal  axis 
polarizabilities. 
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Using  the  dipole  response  model  we  can  calculate  the  polarizability  from  EMI  data  collected 
over  a  target  [1].  The  target  is  classified  as  munitions  (or  other  target  of  interest)  or  clutter 
depending  on  whether  its  polarizability  is  more  munitions-like  or  clutter-like.  Conventional 
dipole  inversion  uses  an  iterative  search  procedure  to  determine  the  dipole  response  model 
parameters  which  produce  the  closest  match  between  the  model  and  the  measured  response.  The 
process  is  illustrated  as  a  flow  diagram  in  Figure  1  below. 


Figure  1.  Flow  diagram  for  conventional  dipole  inversion. 

This  procedure  was  originally  developed  for  data  collected  at  various  locations  over  a  target  with 
a  monostatic  (a  single  pair  of  co-located  transmit  and  receive  coils)  EMI  sensor  such  as  the 
Geonics  EM6 1 .  The  new  generation  of  EMI  sensors  purpose-built  for  classification  are  fixed 
arrays  of  transmit  and  receive  coils  which  excite  and  observe  the  target  over  a  broad  range  of 
angles.  This  is  illustrated  in  two  dimensions  for  a  pair  of  different  coil  configurations  in  Figure  2. 
The  green  and  blue  lines  are  field  lines  for  a  pair  of  coils  A  and  B.  The  pairs  of  arrows  at  the 
target  show  the  directions  of  the  primary  field  vectors  Ha  and  Hb  from  the  two  coils. 


Figure  2.  Multi-coil  arrays  exciting  target  with  different  primary  field  directions. 

In  general,  the  fields  from  three  coils  arranged  in  some  reasonable  fashion  can  be  combined  to 
create  field  components  of  unit  strength  in  three  orthogonal  directions  at  the  target  location.  The 
coefficients  of  the  linear  combinations  are  solutions  of  the  linear  vector  equation 
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a(x)HA  +  b(x)HB  +  c(x)Hc  =  I 


where  I  is  the  identity  matrix  and  the  symbol  (x)  signifies  the  outer  product  of  the  coefficient  and 
field  vectors.  This  is  actually  three  sets  of  equations  wherein  linear  combinations  of  fields  H 
from  coils  A,  B  and  C  with  coefficients  a,  b  and  c  produce  unit  vectors  nx,  nx  and  nx  in  x,  y  and  z 
directions  at  the  target,  e.g. 


^x^Ax  T  bxHBx  +  cxHcx  1 

^x^Ay  +  t>xHBy  +  ^xHcy  —  0 

^x^Az  +  bxHBz  +  cxHCz  —  0 

for  the  x  direction.  A  well-conditioned  H  field  component  matrix 


•Ax 

HBx 

Hex 

Ay 

HBy 

Hey 

^Az 

Hcz 

Hcz 

implies  good  coil/position  combination  for  target  illumination  or  observation. 

The  signal  from  a  single  transmit/receive  coil  pair  is  S  =  R  •  BT.  Given  a  set  of  three 
transmit  and  three  receive  coils,  we  can  solve  for  coefficients  pij  and  qij  for  transmit  and  receive 
coil  combinations  which  synthesize  orthononnal  transmit  field  and  receive  pseudo-field 
components  at  the  target  location  and  then  directly  calculate  the  components  of  the  polarizability 
tensor 


By  =  ^PikS(TxkfRx1)qjl. 

k,l 

This  presents  an  alternative  to  conventional  dipole  inversion.  So  long  as  we  use  the  correct  target 
location  we  should  get  the  same  polarizability  components  for  any  transmit/receive  combination. 
The  idea  then  is  to  search  out  the  target  location  which  minimizes  the  dispersion  of  polarizability 
estimates  from  different  sets  of  transmit  and  receive  coils.  Significantly,  the  minimum  dispersion 
provides  a  direct  measure  of  the  uncertainty  in  our  polarizability  estimate. 

Minimizing  Polarizability  Dispersion 

As  an  expedient,  we  use  the  primary  invariant  of  the  polarizability  tensor  for  our  calculations. 
There  are  three  rotationally  invariant  quantities  associated  with  the  polarizability  tensor  [2].  The 
trace  Tr(B)  is  equal  to  the  sum  of  the  diagonal  elements,  and  because  it  is  a  rotational  invariant, 
also  equal  to  the  sum  of  the  eigenvalues.  The  other  two  invariants  are  the  determinant  of  B  and 
the  sum  of  the  detenninants  of  the  2x2  diagonal  sub-matrices  of  B.  By  virtue  of  their  rotational 
invariance,  these  are  equal  respectively  to  the  product  of  the  eigenvalues  and  the  sum  of  products 
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of  the  eigenvalues  taken  pairwise.  The  invariants  are  the  coefficients  of  the  characteristic 
polynomial  of  B,  and  the  principal  axis  polarizabilities  can  be  directly  calculated  from  them  if 
desired.  An  average  polarizability  <(1>  can  be  calculated  by  averaging  the  polarizability  tensor  B 
over  all  orientation  angles.  The  average  polarizability  is  a  scalar  function  of  time  that  describes 
the  overall  average  response  of  the  target.  It  is  equal  to  the  average  of  the  principal  axis 
polarizabilities,  and  hence  equal  to  'A  times  the  trace  of  B  [3]. 

The  basic  algorithm  is  illustrated  schematically  in  Figure  3.  We  iteratively  search  out  the  target 
location  which  minimizes  the  dispersion  in  estimates  of  <f:S>  for  a  variety  of  coil  combinations. 
Once  the  target  location  has  been  determined  we  can  go  on  to  calculate  the  other  invariants  or  the 
full  polarizability  tensor  if  desired. 


Figure  3.  Flow  chart  for  alternative  approach  to  estimating  the  polarizability  of  a  target. 

There  are  N!/[3!  x  (N-3)!]  ways  of  combining  N  transmit  or  receive  coils  in  sets  of  three.  The 
2x2  TEM  array  [4]  has  four  transmit  coils  which  yield  four  Tx  sets  and  four  three-axis  receiver 
cubes  (12  coils  in  all)  which  yield  220  Rx  sets.  Thus  with  this  array  we  could  calculate  the 
polarizability  880  different  ways.  Not  all  of  the  possible  ways  are  equally  viable,  and  there  is 
clearly  redundancy  since  the  2x2  array  produces  only  48  independent  measurements  of  the  EMI 
response. 

Our  processing  algorithm  selects  a  subset  based  on  the  signal  strength  and  how  well  the  receivers 
interrogate  the  nominal  target  location.  These  are  determined  by  the  norms  and  the  condition 
numbers  of  the  field  component  matrices  of  the  different  possible  combinations.  The  norm  |  |H|  | 
of  the  field  component  matrix  H  is  a  gauge  of  the  field  strength  at  the  target.  We  use  the 
Euclidian  (L2)  norm,  given  by  the  largest  eigenvalue  of  H.  It  is  proportional  to  the  loop  area  and 
the  number  of  turns,  and  more  or  less  inversely  proportional  to  the  third  power  of  the  range  to  the 
target.  The  condition  number  k—  |  |H|  |  |  |H_1|  |,  which  is  equal  to  the  ratio  of  the  largest  to  the 
smallest  eigrnvalue,  is  a  gauge  of  coil  combination  quality  for  target  illumination/observation. 
When  k=  1  all  principal  axes  are  being  interrogated  equally  well.  Large  values  of  Vindicate  that 
at  least  one  of  the  axes  is  not  being  interrogated  very  well.  Larger  |  |H|  |  and  smaller  jkmake  for 
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better  coil  combinations.  Figure  4  shows  the  distribution  of  field  matrix  norms  and  condition 
numbers  for  the  2x2  receiver  array  with  a  shallow  target  near  the  center  of  the  array. 


Target  x  =  -1 .5  cm,  y  =  3  cm,  d  =  5  cm 


Figure  4.  2x2  receiver  array  field  matrix  norm  and  condition  number  for  a  shallow  target  near  the 
center  of  the  array. 

We  use  a  subset  of  the  880  possible  transmit/receive  sets  (typically  -25-50)  with  the  largest 
norms  and  smallest  condition  numbers  for  the  nominal  target  location  in  our  calculations  -  i.e. 
those  to  the  lower  right  in  Figure  4.  The  error  surface  for  target  location  using  our  approach  is 
appears  to  be  well  behaved.  Figure  5  shows  the  polarizability  dispersion  error  surface  over  the 
x-y  plane  at  the  target  depth  for  target  BE-529  (a  60mm  mortar)  from  the  2011  Camp  Beale 
demonstration  [5].  Here  we  use  one  set  of  three  transmit  coils  and  25  sets  of  receive  coils,  and 
calculate  the  dispersion  of  the  net  polarizability  P  =  Ei.jBii(tj)  as  the  ratio  of  the  standard 
deviation  of  the  P  values  for  the  25  receiver  combinations  to  the  mean,  a(P) /P.  In  this  example, 
we  sum  over  all  time  gates  tj  >  88  fj. s.  The  minimum  of  the  error  surface  provides  a  direct 
measure  of  the  uncertainty  in  the  polarizability  estimate. 


0.30 


Figure  5.  Polarizability  dispersion  error  surface  for  target  location  for  Camp  Beale  target  BE-529. 
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We  use  the  basic  Nelder  and  Mead  downhill  simplex  algorithm  [6]  to  search  for  the  target 
location  which  minimizes  the  polarizability  dispersion.  Once  the  target  location  is  determined  we 
calculate  the  principal  axis  polarizabilities.  Figure  6  compares  principal  axis  polarizabilities 
calculated  by  minimizing  polarizability  dispersion  (lines)  with  polarizabilities  from  standard 
dipole  inversion  (points)  for  two  of  the  Camp  Beale  targets.  BE-407  is  a  horseshoe  and  BE-529 
is  a  60  mm  mortar.  For  BE-529,  the  dipole  inversion  calculation  took  4.02  s  using  IDL  on  a 
Lenovo  Thinkpad  X201,  while  the  polarizability  dispersion  minimization  took  3.12  s.  The  effects 
of  signal-to-noise  ratio  (SNR)  and  background  leveling  on  convergence  will  be  analyzed  using 
data  from  the  Camp  Beale  demonstration. 


Figure  6.  Principal  axis  polarizabilities  calculated  by  minimizing  polarizability  dispersion  (lines) 
compared  with  polarizabilities  from  standard  dipole  inversion  (points) 

Signal  to  Noise  and  Background  Leveling 

With  the  2x2  TEM  array,  noise  and  background  leveling  are  intertwined.  At  early  times  the 
background  response  is  dominated  by  coil  ring-down  effects,  while  at  later  times  electromagnetic 
noise  effects  become  more  important.  Background  (no  target)  signals  are  subtracted  from  data 
collected  over  a  target  before  processing.  When  the  sensor  is  being  used  in  the  static  mode  for 
cued  data  collection,  standard  practice  is  to  take  background  shots  every  20-30  minutes. 
Variability  of  the  background  signals  provides  an  estimate  of  the  noise  level.  We  have  chosen  a 
subset  of  the  2x2  TEM  data  from  the  Camp  Beale  demonstration  for  analysis.  Figure  7  shows  the 
site.  The  boxes  show  areas  where  2x2TEM  data  were  collected,  and  the  dots  show  the  locations 
which  appeared  to  be  target- free  and  were  used  for  background  shots.  Our  analysis  uses  data 
collected  on  June  8,  201 1  in  the  open  area  outlined  in  red  with  green  dots.  There  were  14 
background  shots  spread  out  over  a  6-hour  period.  Figure  8  shows  the  mean  background  level  for 
the  various  transmit/receive  combinations.  The  monostatic  (z  axis  receivers  along  the  block 
diagonal)  combinations  show  the  most  pronounced  ringdown.  Standard  deviations  of  the 
background  signals  for  the  various  transmit/receive  combinations  are  shown  in  Figure  9.  The 
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apparent  noise  for  the  monostatic  combinations  is  a  bit  different  from  the  others,  possibly 
reflecting  some  contribution  from  soil  response  variations  over  the  region. 


647200  647300 _  647400  _  647500  _  647600 


Figure  7.  Ariel  photograph  of  Camp  Beale  site  showing  2x2  TEM  data  collection  areas  and 
background  shot  locations. 

Data  were  collected  over  152  targets  during  this  time  period.  Figure  10  shows  the  distribution  of 
signal  to  noise  ratio  (SNR)  vs.  decay  time  for  these  targets.  The  median  and  interquartile  ranges 
are  shown  by  the  red  curves.  SNR  is  calculated  as  the  average  over  all  of  the  48  signal  channels 
(four  transmitters  and  four  three-axis  receivers)  of  the  background  subtracted  signal  divided  by 
the  corresponding  noise  estimate  from  Figure  9.  Note  that  the  SNR  peaks  at  about  0.1  ms. 

We  tested  convergence  of  the  new  algorithm  on  these  data  using  as  a  metric  the  number  of 
function  calls  to  calculate  the  polarizability  dispersion  from  the  downhill  simplex  algorithm.  The 
results  are  given  in  the  next  section. 


8 


mean  background  (mV/A)  mean  background  (mV/A)  mean  background  (mV/A)  mean  background  (mV/A) 


TxO/RxO 


TxO/Rxl 


Tx0/Rx2 


Tx0/Rx3 


time  (ms) 
Tx2/Rx0 


time  (ms) 


< 

> 

E 


> 

E 


Tx2/Rx1 


0.1  1.0  10.0 


time  (ms) 
Tx3/Rx1 


time  (ms) 


< 

> 

E 


Tx1/Rx2 


time  (ms) 


$ 

> 

E 


Tx2/Rx2 


time  (ms) 


< 

> 

E 


Tx3/Rx2 


time  (ms) 


< 

> 

E 


0.1  1.0  10.0 


time  (ms) 
Tx3/Rx3 


Figure  8.  Mean  background  levels  for  2x2  TEM  data  collected  on  June  8,  2011. 
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Figure  9.  Standard  deviation  of  background  levels  for  2x2  TEM  data  collected  on  June  8,  2011 
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Figure  10.  Signal  to  noise  ratio  vs.  decay  time  for  2x2  TEM  target  data  collected  at  Camp  Beale  on 
June  8,  2011. 
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Results  and  Discussion 


We  tested  convergence  of  the  new  algorithm  using  the  Camp  Beale  data.  We  only  used  those 
anomalies  for  which  only  one  target  was  found  during  the  post-test  intrusive  investigation.  This 
amounted  to  106  of  the  original  152  anomalies.  The  downhill  simplex  algorithm  needs  a  fairly 
good  starting  point,  so  we  tested  using  initial  guesses  for  the  target  location  under  each  of  the 
four  quadrants  of  the  array,  at  x  =  + 1 5  cm  and  y  =  + 1 5  cm  (the  centers  of  the  sensors  in  the 
array  are  at  x  =  +20  cm  and  y  =  +20  cm)  and  at  15  cm  depth.  Polarizability  dispersion  was 
calculated  over  all  four  transmit  combinations  and  twelve  receive  combinations,  averaged  over 
time  gates  ranging  from  0.1 -0.6  ms.  Using  the  four  starts,  the  algorithm  converged  for  all  of  the 
targets.  Most  of  the  anomalies  (66  of  the  106)  converged  from  all  of  the  quadrants,  while  only 
two  converged  from  just  one  quadrant.  Of  the  remaining  anomalies,  31  only  converged  from 
three  of  the  four  quadrants  and  seven  converged  from  only  two.  The  number  of  function  calls  to 
calculate  the  polarizability  dispersion  is  a  gauge  of  how  rapidly  the  algorithm  converges,  and  is 
plotted  as  a  function  of  SNR  for  the  106  single-target  anomalies  in  Figure  11.  All  four  quadrant 
starts  are  included,  with  a  limit  of  150  calls  each.  If  a  target  converged  for  none  of  the  quadrant 
starts  there  would  be  600  calls.  The  median  number  of  calls  per  target  is  382,  and  there  is  no 
obvious  trend  with  SNR. 


Figure  11.  Function  calls  to  calculate  polarizability  vs.  signal  to  noise  ratio. 

The  minimum  polarizability  dispersion  a(P)/P  does  vary  with  SNR  (Figure  12,  left).  SNR  is 
calculated  as  the  average  over  all  of  the  48  signal  channels  of  the  background  subtracted  signal 
voltage  divided  by  the  corresponding  noise  estimate  from  Figure  9.  As  one  might  expect,  targets 
with  higher  SNR  tend  to  have  less  uncertainty  (dispersion)  in  the  polarizability  estimate  than 
those  with  lower  SNR.  Although  there  is  a  good  bit  of  spread,  the  polarizability  dispersion  based 
inversion  data  seems  to  be  reasonably  well  fit  by  a  straight  line  in  log-log  space  of  the  form 
ct(P)/P  oc  SNR-1,  consistent  with  the  Cramer-Rao  bound  on  parameter  estimation  error  [7]. 
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Figure  12.  Residual  error  vs.  signal  to  noise  ratio  for  polarizability  dispersion  and  standard  dipole 
inversion. 


Fit  error  for  conventional  dipole  inversion  (basically  a  {data  —  model) /data  )  generally 
follows  a  similar  trend  [8],  but  there  is  a  bit  more  scatter  in  the  plot  of  dipole  fit  error  vs.  SNR 
for  these  data  (Figure  12,  right)  than  in  the  corresponding  polarizability  dispersion  plot  on  the 
left.  There  is  little  correlation  (r  =  0.22)  between  the  minimum  polarizability  dispersion  and  the 
mismatch  between  data  and  model  with  conventional  dipole  inversion.  This  most  likely  reflects 
the  fact  that  the  two  approaches  treat  the  data  differently,  but  the  reason  for  the  increased  scatter 
with  the  dipole  fit  error  remains  unclear. 


Standard  Dipole  Inversion 


With  conventional  dipole  inversion  the  standard  metric  for  data-model  match  is  the  mean 
squared  deviation,  in  which  case  it  is  a  separable  nonlinear  least  squares  problem  [9].  The  signal 
S  depends  linearly  on  the  polarizability  B,  but  the  transmit  field  T  and  receive  pseudofield  R 
used  to  calculate  the  signal  both  depend  nonlinearly  on  the  target  location  relative  to  the  sensor. 
The  usual  approach  with  such  problems  is  to  iteratively  search  on  the  nonlinear  parameters 
(target  x,  y,  z  location  in  this  case)  while  solving  the  linear  least  squares  problem  for  the  linear 
parameters  (the  elements  of  the  polarizability  tensor)  at  each  step.  With  the  2x2  TEM  array  there 
is  an  overdetennined  set  of  48  linear  equations  in  six  unknowns  (B  is  symmetric  so  only  six  of  its 
nine  elements  are  unique)  to  be  solved  at  each  step.  The  solution  minimizes  the  mean  squared 
error  in  matching  all  of  the  data.  With  polarizability  dispersion  based  inversion  we  use  a  similar 
tactic  of  iteratively  searching  for  the  target  location  which  minimizes  an  objective  function,  but 
now  that  objective  function  reflects  the  uncertainty  in  the  polarizability  estimate  rather  than  the 
actual  mismatch  between  observed  and  modeled  signals.  The  polarizability  calculations  are 
equivalent  to  solving  6-equation  subsets  of  the  48  which  use  three  different  transmitters  and  three 
different  receivers.  We  are  clearly  weighting  the  data  differently  with  the  two  inversion 
approaches,  and  should  expect  somewhat  different  behavior. 

In  at  least  some  cases  we  see  improved  performance  using  the  polarizability  based  objective 
function.  Figure  13  shows  results  of  test  stand  measurements  of  a  37  mm  projectile  with  the  2x2 
TEM  array.  The  plot  shows  the  errors  in  calculated  polarizabilities  for  the  two  processing 
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techniques  for  increasing  vertical  separation  between  the  sensors  and  the  projectile.  The 
polarizability  error  plotted  in  this  figure  is  the  absolute  difference  between  the  calculated  (3i(t) 
and  reference  37  mm  (3;(t)  at  62  time  gates  ranging  from  0.09  ms  to  2.77  ms,  normalized  by  the 
reference  polarizability.  For  the  reference  (3s  we  use  the  average  of  the  standard  dipole  inversion 
(3s  and  the  polarizability  based  (3s  at  29  cm.  Beyond  50  cm  we  see  an  obvious  increase  in  the  (3 
error  for  standard  dipole  inversion  (black  diamonds)  relative  to  our  polarizability  based 
processing  (blue  circles),  although  both  are  increasing  rather  rapidly  at  this  point.  The 
improvement  is  due  to  the  fact  that  the  polarizability  based  processing  is  getting  a  better  location 
estimate.  But  the  advantage  is  soon  lost.  At  greater  depths  the  target  is  simply  not  being 
illuminated  well  enough  by  the  sensor. 


Distance  below  array  (m) 

Figure  13.  Errors  in  calculating  the  polarizability  of  a  37  mm  projectile  at  various  distances  below 
the  2x2  TEM  array  with  standard  dipole  inversion  and  polarizability  dispersion  based  processing. 
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Conclusions  to  Date 


The  intent  of  this  research  project  is  to  explore  alternatives  to  conventional  dipole  inversion  for 
extracting  target  features  from  multi-axis  EMI  sensor  array  data.  Conventional  dipole  inversion 
searches  for  the  target  parameter  values  (location  and  polarizabilities)  which  minimize  the 
difference  between  measured  signals  and  those  calculated  using  the  dipole  response  model.  Here, 
we  have  considered  an  alternative  approach  that  seeks  to  determine  the  parameter  values  which 
minimize  an  objective  function  based  on  the  dispersion  in  estimates  of  the  target’s  polarizability 
using  different  combinations  of  transmitters  and  receivers.  A  fortiori  the  approach  allows  for 
direct  calculation  of  the  uncertainty  in  the  estimated  polarizability. 

This  interim  report  documents  results  on  the  convergence  properties  of  a  downhill  simplex  based 
algorithm  for  detennining  a  target’s  location  (and  polarizabilities)  using  this  approach. 
Specifically  we  examined  convergence  of  the  algorithm  for  data  collected  at  the  former  Camp 
Beale  in  201 1  and  found  no  impact  due  to  signal-to-noise  ratio  (SNR)  and  background  leveling 
effects.  The  algorithm  converges  properly  for  all  field  targets  irrespective  SNR.  However,  the 
minimum  polarizability  dispersion  does  vary  systematically  with  SNR  -  targets  with  higher  SNR 
tend  to  have  less  uncertainty  (dispersion)  in  the  polarizability  estimate  than  those  with  lower 
SNR. 

Testing  shows  that  polarizability  dispersion  based  inversion  can  produce  more  accurate 
polarizabilities  than  conventional  dipole  inversion  in  some  cases.  However,  the  fraction  of 
targets  in  the  ESTCP  classification  demonstrations  that  cannot  be  accurately  analyzed  using 
conventional  dipole  inversion  may  be  small  enough  (<5%)  that  modest  improvements  in 
calculating  the  polarizability  are  likely  to  have  very  little  effect  on  classification  perfonnance.  In 
the  next  phase  of  this  project  we  will  be  exploring  improved  procedures  for  target  classification 
based  on  principal  axis  polarizabilities. 
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